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Abstract 

We study the chemical nature of the bonding of an oxide layer to the parent metal. In 
order to disentangle chemical effects from strain/misfit, Ti(1010)/TiO2(100) interface has been 
chosen. We use the density functional pseudopotential method which gives good agreement 
with experiment for known properties of bulk and surface Ti and T1O2. Two geometries, a 
film-like model (with free surface in the structure) and a bulk-like model (with no free surface 
in the structure) are used to simulate the interface, in each case with different terminations 
of Ti and Ti02- For the single-oxygen interfaces, the interface energy obtained using these 
two models agree with each other; however for the double-oxygen ones, the relative stability 
is quite different. The disturbance to the electronic structure is confined within a few atomic 
layers of the interface. The interfacial bonding is mainly ionic, and surprisingly there is more 
charge transfer from Ti to O in the interface than in the bulk. In consequence the Ti/Ti02 
interface has stronger binding than the bulk of either material. This helps to explain why the 
oxide forms a stable, protective layer on Ti and Ti alloys. 

INTRODUCTION 

In contrast to the many atomistic studies of metals and oxides, both at the ab initio and empirical 
level, the metal-oxide interface is much less well studied. This may appear surprising given its 
central role in corrosion. One reason may be the difficulty in identifying a candidate interface 
where the chemical effects are not entangled with misfit strain energy. Another is the relatively 
large size of system which is required to isolate the system. Although the strain energy is long- 
ranged, there are reasons to expect that the chemical effects of the interface may not be: on the 
metal side the free electrons should screen the Coulomb forces, while the image charges induced in 
the metal should mimic the electrostatics of a bulk oxide. At the interface itself, such classical ideas 
break down, and a full quantum treatment is needed to determine the nature of the bonding. In this 
paper we will consider the low-misfit Ti(1010)/TiO2(100) interface, with a view to determining the 
range over which chemical effects are significant and the nature of the bonding cross the interface. 
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Titanium alloys are widely used in many fields such as the aerospace industry, chemical plants, 
and even sporting goods. The high strength-to-weight ratio and excellent corrosion resistance 
account for its wide application. The corrosion property is mainly a result of the formation of stable 
protective oxide film, which consists primarily of Ti02- However, above 600°C, the fast diffusion 
of oxygen through the oxide layer into the bulk can result in excessive growth of oxide layer and 
embrittlement of the adjacent oxygen rich layer of the titanium alloy, which limits its maximum use 
temperature.- Alloying and coating have been found effective to address this problem. Obviously, 
the interface between the titanium and its oxide plays a vital role in the corrosion, both through its 
adhesive strength and the diffusion of species (O and/or Ti) through it, and its structure is a key 
aspect for understanding the behavior of titanium alloys. 

Extensive experiments on pure titanium 2-5 and titanium alloys 6,7 have established that the 
crystal structure of the oxide is normally rutile (tetragonal, P42/mnm). Although Guleryuz et al.— 
reported some diffraction angles consistent with the anatase structure(tetragonal, I4i/amd), in the 
scale of Ti-6A1-4V oxidized at 600°C, rutile is the dominant oxide at 650°C. Due to the competition 
between surface free energy and strain energy, the growth of rutile on pure titanium exhibits a 
preferential direction with a specific crystallographic orientation relationship (COR) between 
titanium and rutile. Three possible CORs between Ti and rutile were established by Flower et 
al.- 2 - using an in situ method: Ti(0001)[l 120] // TiO 2 (010)[001], Ti(10T0)[0001] // TiO 2 (100)[010] 
and Ti(1120)[0001] // TiO 2 (001)[100]. Among these three different matchings, Ti(10l0)[0001] // 
TiO2(100)[010] has the smallest mismatch on the plane of the interface, which means only a tiny 
interface strain would be required. 

The structure of the oxide and the orientation relationship can be easily determined by the ex- 
periments; however, some other quantities like the chemical composition, atomic structure of the 
interface, and the nature of the bonding (ionic/metallic/covalent) across the interface are currently 
experimentally inaccessible. Fortunately, nowadays, it is possible to deal with coherent interface 
structures (large system, low symmetry) using accurate first-principles theoretical methods to ob- 
tain those quantities. Coherency implies that one part (metal or oxide) will be strained to match the 
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other one perfectly to maintain the coherency without misfit dislocations. In this case, DFT super- 
cell method can be a good tool to study interfaces with small mismatch, but even when the misfit 
is quite big, it has been assumed that the interfacial regions between the misfit dislocation being 
modelled. 9 Such a first-principles supercell method has been used to study many different inter- 
faces, for example, interface of Zr0 2 /Ni,^ AVA\ 2 3 ^ Ti/TiN,— YSZ/AI2O3 ,— Nb/Nb 5 Si 3 ,-^ etc. 
Here we present a DFT study of Ti(10T0)[0001] // TiO 2 (100)[010] interface. The purpose of this 
study is to determine the optimal atomic structure and energy of Ti(1010)[0001] // TiO2(100)[010] 
interface, and characterize the nature of the interfacial bonding. 

METHODS 

The Vienna Ab initio Simulation Package (VASP)— utilizing a plane-wave basis set for the 
expansion of the single-particle Kohn-Sham wave functions, was used in this study. The projector 
augmented wave (PAW) method,— was employed to describe the electron-core interaction. The 
3p semicore electrons of Ti were treated as valence, given a 10-electron PAW-pseudopotential. 
For the exchange-correlation interaction, we adopted generalized gradient approximation (GGA) 
as parameterized by Perdew and Wang (PW91).— A high cutoff energy of 525 eV was used. 
Sampling of the Brillouin zone was performed with a Monkhorst-Pack grid.— Ground-state atomic 
structures were obtained by minimizing the Hellmann-Feynman forces on the atoms, and all the 
atoms were free to relax. The relaxations terminate when the maximum force on the atoms is less 
than 0.05 eV/A. For some calculations a dipole moment is present, in such cases the divergent 
terms are removed by the Ewald sum; we did not include a dipole correction 7 as previous work in 
TiC>2 has shown this to have minor effects. 
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Table 1: Structure parameters (internal coordinate u in TiCh refers to the position of the 
oxygen atom in the unit cell), bulk modulus Bq (GPa), compared with experimental data. 
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Figure 1 : DOS per TiC>2 unit in bulk rutile. 



RESULTS: parameters and basic properties of Ti and TiC>2 
Bulk properties 



To verify the accuracy of our computational parameters, we first calculated the bulk properties of 
Ti and rutile. The k-point mesh was set at 5 x 5 x 8 and 1 1 x 1 1 x 7 for the bulk TiC>2 and hcp-Ti, 
respectively. These provide convergence to within 1 meV and as can be seen from \T\ the calculated 
lattice parameter and bulk modulus data agree excellently with the experiments. From the DOS of 
Ti02 (©> the gap between highest occupied and lowest unoccupied state of Ti02 is about 1 .7 eV, 
far away from the experimental value, which is around 3.0 eV:— this does not affect the bonding 
and is typical of the error of DFT in describing excited states. 

Surface properties 

To determine the minimum slab thickness needed to reliably calculate the surface/interface we ex- 
amined the convergence of the relaxation and energy of the surfaces with respect to slab thickness. 
In this calculation, we used k-point sampling 8x5x1 and 5 x 8 x 1 for Ti(1010) and TiO2(100). 
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In both slabs, a 15 A vacuum region is introduced to avoid the interaction between periodic im- 
ages. For both surfaces, only the (1 x 1) structures (along the two lattice vectors of the surface, the 
symmetry is the same as in the bulk) are considered, since these have been observed by for clean 
surface under normal conditions by experimentalists.— - — 

The calculated surface energy might diverge with the thickness of the slab, if there is any 
numerical difference between the calculation of bulk and slab, arising, for example, from the k- 
point mesh. Following Boettger— we avoid this possibility, by evaluating the surface energy 7 as 
follows: 



7=( J E s I ? ab -NA J Ei/AN)/2S 

(1) 

A pN' _ pN' _ pN'-AN 
' LVC slab — ^slab ^slab 

where N' is the minimum number of the slab layers for which the energy converges, £^ ab is the 
energy of a N-layer slab, and S is the area of the surface. 

Ti(lOlO) 

In hep Ti, two different (1010) planes exist, depending on whether the surface terminates in a 
large (1.694 A) or small (0.847 A) interlayer spacing. From the experimental results, 2 ^ Ti(1010) 
surface with a small first interlayer spacing was the favored one. Our calculation confirmed this, 
and henceforth, we concentrate on the more stable Ti(1010) surface. 

Table 2: Interlayer distance relaxation of the Ti(1010) surface vs slab thickness, shown as a 
percentage of the bulk spacing. The surface energy 7 obtained with different slab thickness 
is also listed, in J/m 2 . 

Layer Mn AJ23 AJ34 AJ45 AJ56 AJ78 AJs9 7 

~l2 -3.26 -5.43 5.98 -2.20 3.65 L99~ 

14 -3.64 -4.74 4.10 -1.13 2.19 2.16 1.99 

16 -3.75 -4.80 3.95 -1.01 1.97 1.25 -0.29 1.99 

18 -3.70 -4.93 3.99 -1.18 2.07 1.39 -0.39 1.99 

Expt. fl -5.9 



ref24 
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The surface relaxation data is shown in [21 the distance between the first two layers d\2 de- 
creases by about 4% compared with that in the bulk, in reasonable agreement with 6% from the 
experiment— (the discrepancy is less than 0.02 A). J23 contracts and J34 expands, showing the 
same trend as the results by self-consistent tight bonding model. 27 Relaxations between deeper 
layers are even smaller (less than 2%). The surface energy converged quickly with cell size to 1.99 
J/m 2 . This is slightly higher than the 1.91 J/m 2 of Ti(0001) surface which we obtained using the 
same method: Ti(0001) has an experimental energy of 1.99 J/m 2 . 28 Thus Ti(1010) is just slightly 
less stable than the close-packed Ti(0001) surface. For further study, we chose a 16-layer slab to 
simulate Ti(10T0). 

TiO 2 (100) 

Extensive theoretical studies have been done to investigate the TiC>2(100) surfaces. There are three 
possible terminations for TiO2(100), which we call O-Ti-O, O-O-Ti, and Ti-O-0 (named by the 
atomic arrangement from the surface to inner layer). The first of these © has the smallest surface 
polarization due to the symmetrical arrangement of O, and the geometry has been observed by 
experiment.— We find it to be the most stable one, and we use the O-Ti-0 terminated TiO2(100) 
surface to determine the minimum thickness of the slab needed. One thing that should be noted is 
that in this subsection, the O-Ti-0 unit is treated as one layer, i.e., the 5-layer slab mentioned here 
has 15 atomic layers. 

The choice of exchange-correlation functional has been found to affect the surface energy quite 
significantly, 29 and normally LDA gives a higher surface energy than GGA. Due to the lack of 
the experimental results, it is difficult to clarify which function best describes the real situation. 
Here, we compared our results with other calculations using the same exchange-correlation func- 
tion (GGA), in order to verify the accuracy of our calculation. As shown in [31 Oi, T12, and O3 
exhibit large relaxations along [010]. Moreover, Ti2 relaxes inward along [100], while outward 
relaxation is observed for Oi and O3. These relaxations increase the effective coordination of Ti2 
(fivefold).—^ Our surface relaxation results agree with that found by Muscat et al.. 29 14] lists the 
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surface energy dependence on the slab thickness: it converges quickly to 0.68 J/m 2 which agrees 
very well with Perron et al.— (10 valence electrons are considered for Ti) and Labat et al. (GGA- 
PBE), 30 while it is smaller than Muscat's GGA result. 2 - Considering the convergence of both the 
surface structure and energy, we believe that a 13-layer slab is thick enough to model TiO2(100). 



V 



V 



[100] 



[010] 



[001] 



Figure 2: TiO2(100) surface structure. Red spheres represent O atoms, and silver spheres represent 
Ti atoms. 
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Table 3: The displacement of ions (A) (as labeled in|2]> relative to the bulk terminated struc- 
ture of TiC>2(100) surface obtained with different slab thickness. 
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Ref. fl 



Label 


[010] 


[100] 


[010] 


[100] 


[010] 


[100] 


[010] 


[100] 


[010] 


[100] 


1 


0.32 


0.07 


0.32 


0.08 


0.28 


0.06 


0.28 


0.06 


0.31 


0.06 


2 


-0.10 


-0.03 


-0.11 


-0.02 


-0.15 


-0.03 


-0.15 


-0.03 


-0.11 


-0.04 


3 


0.17 


0.06 


0.18 


0.08 


0.14 


0.06 


0.13 


0.05 


0.13 


0.02 


4 


0.07 


0.01 


0.08 


0.02 


0.04 


0.01 


0.04 


0.00 


0.03 


0.00 


5 


-0.05 


0.02 


-0.06 


0.04 


-0.10 


0.03 


-0.10 


0.02 






6 


0.01 


0.03 


0.04 


0.04 


0.00 


0.02 


0.00 


0.02 






7 


0.01 


-0.01 


0.04 


0.01 


0.01 


-0.01 


0.01 


-0.01 






ref29. 























Table 4: The dependence of surface energy of TiC>2(100) on slab thickness, in unit of J/m 2 . 

Layer 5 9 13 17 Ref° Ref* Ref. c 



reflail; b ref|30|; c ref|29 
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0.70 0.68 0.68 0.68 0.67 0.69 0.83 
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RESULTS: Interface properties 

The orientation relationship is set as Ti(1010)[0001] // TiO2(100)[010] across the interface. The 
size of the Ti(lOlO) surface cell is 2.934 A x 4.638 A, while it is 2.971 A x 4.645 A for 
TiO2(100). A coherent interface is obtained by a small strain of the softer Ti to perfectly match the 
TiO2(100) surface cell: this has little effect on the energy since the mismatch of the two surface 
cells is so small. To check dependence on boundary conditions we used two models to simulate 
the interface, a film-like model with an interface and two free surfaces, and a bulk- like model with 
two interfaces. The interface structures were relaxed with k-point mesh 3x3x1 (four irreducible 
k points). 

Film-like interface model 

The film-like model is generated with 16-layers of Ti, 13-layers of TiC>2 and a 15 A vacuum 
region. This gives one interface and two free surfaces. 

Various terminations of the surfaces and metal/oxide interface are possible, their stability de- 
pending on the environmental condition like the partial pressure of O2 gas.— As we mentioned 
above, TiO2(100) has three possible terminations, with uppermost layers O-O-Ti, O-Ti-O, and 
Ti-O-O. For interfaces, we considered different stacking sequences of Ti onto these terminations. 
Following work on other interfaces,-^ the first layer of the Ti-metal is placed... 

• ' OT' : directly above uppermost Ti cations; 

• ' HCP' : above the second layer of Ti cations; 

• ' FCC : above the third layer of Ti cations. 

• ' TT' : for O-O-Ti, as an extension of the oxide. 

• ' TT' : for Ti-O-O, as an extension of the metal. 

The ' TT' configurations are those which would allow growth of the oxide/metal by simple 
extension of a stable interface, with the Ti ' atoms' and Ti ' cations' directly adjacent. In fact, the 
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two ' TT' configurations are the same at the interface, the calculation differing in the accompany- 
ing surface. 

Despite this apparently exhaustive survey of stacking sequences, we also found that where 
the TiO2(100) surface is quite corrugated, significant reconstruction can occur within the Ti-metal 
region, spontaneously generating a stacking fault. When this happens, we report the most stable 
relaxed structure, and label it ' HCP-2' and ' FCC-2' in [5] 

To avoid spurious contraction due to surface tension, we fix the lattice parameters in the in- 
terface plane. But we allow relaxation perpendicular by first calculating the total energy of the 
unrelaxed interface structure as a function of interface separations d. Once the optimal value is 
found, all the atoms are relaxed to the ground state at fixed volume. 

In order to evaluate the strength of the interface, we calculated the work of adhesion, W a d, 
which is defined as: 

Wad = (Em + E o l ~ e m/o)/ s 

where is the surface energy of Ti-metal, Eq 1 is the surface energy of oxide, E^ Q is the 
energy of the interface structure, and S is the interface area. A positive value of W aa - means that 
the interface is energetically favorable over the free surfaces. [5] shows that W a d increases with 
increasing number of interfacial O atoms (Ti-O-O, O-Ti-O, and O-O-Ti). However, much of this 
is due to differences in surface stability; e.g. the O-O-Ti and Ti-O-0 with ' TT' stacking have the 
same chemical composition near the interface but very different Wad- 
To understand the growth of an additional Ti02 layer, each of these interface types needs to be 
considered. In the following subsections, the interface with the most stable stacking sequence for 
each termination is analyzed in detail. 

For all configurations listed in [61 for the metal slab, the Ti atom in the center has a charge 0.00, 
as in bulk Ti, which again confirms that the slab we used is thick enough. For the Ti02 slab, the 
charge associated with central Ti and O ions are 2.18e and -l.lOe, respectively. This agrees well 
with DFT+U results from Jess and the co-workers— who found the Bader charge is 2.22e and 
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Table 5: Relaxed values of W ac i (J/m 2 ) for interfaces with different termination and stacking 

o 

sequence, with interfacial distance d (A). 









O-O-Ti 






O-Ti-O 






Ti-O-0 






OT 


HCP 


FCC TT 


HCP-2 


OT 


HCP FCC 


FCC-2 


OT 


HCP FCC 


TT 


d 


1.65 


0.57 


1.09 0.97 


1.17 


1.87 


1.37 0.33 


0.31 


1.90 


1.66 1.58 


0.55 


W ad 


7.54 


9.43 


8.99 9.32 


10.34 


2.28 


2.57 3.49 


3.55 


1.41 


2.59 2.04 


2.83 



-X.Yle for Ti and O. Although Bader charge values are far from the formal charge Ae (Ti 4+ ) and 
-2e (0 2 ~), we believe that useful bonding information can be obtained from the charge transfer. 

Atomic, electronic structure of O-O-Ti interface 



o 



• o 

. Ti 4 

Ti, ° 

o o 

. Ti 2 

Tii ° 





TL 



Oi 



Ti 3 Ti 4 



coo 



Ti 5 1 



Figure 3: ' HCP-2' structure of O-O-Ti interface (only some layers near the interface is shown), 
viewed along [001] direction of Ti02, the dashed line indicates the position of the interface. Left: 
without atomic relaxation (the equivalent atoms on the boundary are also shown). Right: with 
atomic relaxation. Large red spheres are O atoms, and small silver spheres are Ti atoms. 



The stacking sequence ' HCP-2' gives rise to the most stable O-O-Ti interface ©. The struc- 
ture is displayed in [3l To observe the geometry change due to the atomic relaxation clearly, the 
unrelaxed (input) structure is also shown and the ions with large displacement are labeled. On the 
metallic side, Ti2 and Ti4 move towards the oxide. Tii and Ti2 are found shifting along [0001] 
of Ti significantly. As a result, the surface morphology of Ti(1010) is totally changed: the Tii 
and Ti2 layers merge to give a flat surface; the spacing between Ti2 and Ti3 layer changes to 2.19 
A from 1 .69 A in the bulk; the interlayer spacing between Ti3 and Ti4 decreases to 0.43 A from 
the bulk spacing 0.85 A. For the oxide slab, the relaxation is less significant, the interfacial atoms 
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move towards metal slab, in particularly O2, thus the surface of the oxide becomes flatter. Across 
the interface, the separations of Tii-Oi and Ti 2 -0 2 is 2.00 A and 2.15 A, which is close to the 
Ti-0 bond lengths in bulk Ti0 2 (1.96 A and 2.00 A). 2 is located in the center of the ' half 
octahedron (formed by Tii, Ti 2 and Tis), and Oi sits in a similar octahedral site. 

Bader charge analysis 34 is applied to study the charge transfer. In[6l we give the Bader charge 
associated with the labeled atoms © as well as the atoms in the center of the metal/oxide slabs, 
and the value of the charge is the variation from the neutral Ti/O atom (lOe for Ti, and 6e for O). 
Thus a negative value means accepting negatively charged electrons, while a positive value means 
donating electrons. 

Table 6: Bader charge (e) of the selected atoms in the various interface structures (denoted 
by the termination, the model used, and the staking sequence in the bracket). Definitions of 
the atoms are given in the appropriate figures. Bulk values for Bader charge are 0.00 in Ti, 
and 2.18e and -l.lOe for Ti and O in bulk rutile. For all the structures listed in the table, the 
Ti atom/ion in the center of the slab has the same Bader charge as in the bulk. 



Ti slab Ti0 2 slab 



Atom 


Tii 


Ti 2 


Ti 3 


Ti 4 


Ti 5 


Oi 


o 2 


O(center) 


O-O-Ti film(HCP-2) 


0.81 


0.82 


-0.21 


0.32 


2.08 


-1.20 


-1.28 


-1.11 


O-Ti-0 film(FCC-2) 


0.82 


0.06 


0.24 


0.01 


1.78 


-1.37 


-1.24 


-1.10 


Ti-O-0 film(TT) 


0.56 


-0.06 


0.18 


-0.04 


1.12 


-1.27 


-1.20 


-1.09 


S-O fl bulk(FCC) 


0.77 


0.15 


0.23 


-0.09 


1.78 


-1.39 


-1.24 


-1.11 


D-O b symmetry bulk(TT) 


1.05 


0.62 


-0.01 


0.12 


2.10 


-1.21 


-1.27 


-1.11 


D-O b asymmetry bulk(TT) 


1.06 


0.60 


-0.05 


0.15 


2.11 


-1.27 


-1.21 


-1.11 


D-O b asymmetry bulk(HCP-2) 


0.83 


0.81 


-0.19 


0.31 


2.08 


-1.19 


-1.28 


-1.11 



a Single-Oxygen; b Double-oxygen. 

For ' HCP-2' each of the interfacial atoms Tii an d Ti 2 donates 0.8 electrons, which account for 
the net charge transfer of 1.61 electrons from the metal slab to the oxide ©. Bader analysis gives 
a negative charge for Ti 3 , but this appears to be an artifact due to the large relaxation around it, 
leading to a large Bader volume, which in turn encloses more electrons. It is better to consider the 
layer containing Ti 3 and Ti4, which has overall a small positive charge. The interfacial oxygen ions 
Oi and 2 attract more electronic charge than the oxygen in the center, due to the more electron 
donators nearby. Tis transfers fewer electrons than the Ti ions in the center of the oxide, which is 
also understandable as its coordination is smaller. The large electron transfer found here implies 
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strong ionic bonding across the interface. 

The electronic density of states is projected onto selected atoms to determine the bonding 
character ©. Notably, Ti3 and T14 are dissimilar, as is reflected in the peak/valley just below the 
Fermi level. The interfacial atoms Tii and Ti2 show a small hybridization peak in the region from 
-7.5 eV to -2.5 eV, representing weak covalent bonding to the oxygen nearby. Strong hybridization 
between the Ti 3d and O 2p is observed from the DOS of all ions on the TiC>2 side of the interface, 
indicating some covalent bonding in rutile.— T15 has almost no Ti conduction band density: which 
implies that bonding in the O-O-Ti structure is not metallic. The DOS of Oi and O2 are similar to 
the oxygen in the center of oxide slab, with a small shift of both s and p peaks to a lower energy 
level indicating a stronger Madelung field at the interface, which can stabilize the system.— For 
the oxygen on the surface, the s orbital shows no shift, and the width of its p orbital is reduced. 
By comparing the total DOS of the interface slab (Oa)) an d DOS of pure Ti02 ® , we see a small 
peak around -7.5 eV induced by the interface atoms. This will be discussed further later. 
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Figure 4: DOS projected onto selected atoms for the O-O-Ti interface. The labels of the atoms are 
identified in [31 Note that within Ti02 all the states corresponding to the 4s3d band are unoccupied, 
consistent with formal ionic charges of Ti 4+ and 2 ~. 



Atomic, electronic structure of O-Ti-O interface 

As shown in [51 for O-Ti-0 interface systems, the interface with strongest adhesion is ' FCC-2' , 
which has a similar structure to ' FCC .[6] shows the optimized ' FCC-2' interface: Tii and Ti3 
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Energy relative to Fermi Level (eV) 

Figure 5: Total DOS of various interfaces using two models, (a) Total DOS of the O-O-Ti ter- 
minated interface using film-like model, (b) Total DOS of the O-Ti-0 terminated interface using 
film-like model, (c) Total DOS of the Ti-O-0 terminated interface using film-like model. In this in- 
terface slab, a O-O-Ti terminated TiO2(100) free surface is included. This polarized surface might 
account for the abnormal shape of the 2s states of oxygen, (d) Total DOS of the single-oxygen 
interface using bulk-like model, (e) Total DOS of the double-oxygen interface using symmet- 
ric structure of bulk-like model (including two ' TT' stacking interfaces), (f) Total DOS of the 
double-oxygen interface using asymmetric structure of bulk-like model (including one ' TT' and 
one'HCP-2' stacking interface). 



Ti 4 



Tu 



9_ °_ 

Figure 6: ' FCC-2' structure of O-Ti-0 interface (only several layers near the interface is shown), 
viewed along [001] direction of Ti02, the dashed line indicate the position of the interface. Left: 
without atomic relaxation (the equivalent atoms on the boundary are also shown). Right: with 
atomic relaxation. Red spheres are O atoms, and silver spheres are Ti atoms. 
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are shifted significantly toward the oxide, forming planes containing both Ti and O atoms. With 
this displacement, the interlayer spacing of the metal slab near the interface changed to 1 . 19 A and 
1.36 A, which are 0.85 A and 1.69 A in the bulk Ti, respectively. Relaxation of Oi and O2 leads 
to a displacement towards the metal even larger than that in the free surface structure. Finally, the 
separation O1-T12 and 02-Tii is 2.13 A and 2.10 A, close to the spacing in bulk Ti02 (1.96/2.00 
A). 

The Bader charge analysis results are listed in [6] The most significant charge transfer comes 
from Tii and Ti3, which donate a total of 1.16 electrons from the metal to the oxide. The excess 
charge is located on the interfacial atoms (Tis, Oi and O2). 

As seen from H the atom-projected DOS converges rapidly to bulklike values away from the 
interface. For the interfacial ion Tis, states in the band gap of the oxide implies a kind of metallic 
behavior, and it is a little stronger than that in O-O-Ti interface, because of its direct exposure to 
the metal slab. The shift of DOS is very apparent in the interfacial atoms Oi and O2, especially for 
Oi, the 2s states is below -20 eV. For the surface oxygen, the 2s and 2p orbitals move to a higher 
energy level, which contribute to the small peak near -17.5 eV in the total DOS Qb)). Also, 
comparing Q] (pure Ti02),Ob) (interface with surface), and[5£d) (interface with no surface), it can 
be concluded that, in[5tb), the small peak below -20 eV and peak near -7.5 eV are associated with 
the interface atoms, while the peak near -17.5 eV comes from surface atoms. These agree with the 
our calculation results that the formation of the interface is exothermic, while that is endothermic 
for the surface. 

Atomic, electronic structure of Ti-O-O interface 

For the Ti-O-0 structures, ' TT' stacking exhibits the lowest energy. The input and relaxed struc- 
ture of 'TT' stacking are shown in[8l Again, the relaxation tends to flatten the surfaces. The 

o o 

interfacial spacing (distance between Tii layer and Ti$ layer) is 0.55 A, smaller than the 0.85 A of 
the interlayer spacing in bulk Ti. Oi moves towards the interface, which increases the effective 
coordination of Ti$. The downward shift of Tii an d Ti3 are also observed. The interlayer spacing 
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Figure 7: DOS projected onto selected atoms for the O-Ti-0 interface. The label of the atom is 
identified in[6l 
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Figure 8: ' TT' structure of Ti-O-0 interface (only several layers near the interface is shown), 
viewed along [001] direction of TiC>2, the dashed line indicate the position of the interface. Left: 
without atomic relaxation (the equivalent atoms on the boundary are also shown). Right: with 
atomic relaxation. Red spheres are O atoms, and silver spheres are Ti atoms. 
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o 

of Tii-Ti2 and Ti2-Ti3 is 1.88 A and 0.74 A, respectively, which slightly deviate from the spacing 
in the bulk Ti. 

The Bader charge analysis © and DOS projection © are qualitatively similar to the two 
systems. However, only OAle is moved to the oxide slab from the metal, much less than in the 
other two surfaces. From the projected DOS of the interfacial ions (Tii,Ti2, Ti3 and T14), no 
covalent features are evident between the two slabs, but metallic bonding can be inferred from the 
occupied 4s3d states in Tis which is adjacent to the metal slab. The relatively small W a d of the 
Ti-O-0 systems can be understood by the pDOS ©, where the 2s states of surface oxygen move 
to a very high level just below -15 eV. Thus the states between -15 eV and -17.5 eV in the[5tc) 
come from the oxygen near the free surface (O-O-Ti terminated) included in the slab model. 
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Figure 9: DOS projected onto selected atoms for the Ti-O-0 interface. The label of the atom is 
identified in[8l 



Bulk-like interface model 

In this model, we consider a periodic arrangement ■ • ■ metal-oxide-metal- • ■ , with no free surfaces. 
One advantage of this model over the previous one is that the effects of surface dipole interactions 
are eliminated. 

Considering the termination of TiO2(100) at the interface, there are two structures with different 
layer ordering 



17 



.... Ti-Ti-Ti-O-Ti-O-O-Ti- • • O-O-Ti-O-Ti-Ti-Ti- • • , 

.... Ti-Ti-Ti-O-O-Ti-O-O-Ti- ■ ■ O-O-Ti-O-O-Ti-Ti-Ti- ■ • . 

The former, denoted as single-oxygen interface, may have either two identical interfaces or differ- 
ent stacking sequences. For the latter, double-oxygen interface, we first calculate a ' TT' configu- 
ration where the first cation (metal) layer is like the extension of the metal (oxide). Subsequently 
we made an 'asymmetric' structure: one interface with 'TT' stacking, the other with ' HCP- 
2' stacking, so that the energy of the double-oxygen interface with 'HCP-2' stacking can be 
obtained. Again, we fix the dimensions of the cell in the interface plane. 

In order to evaluate the strength of the interface, we defined the interface energy Wbulk as: 

where E^ mQn is the energy of interface structure consisting of m Ti atoms and n O atoms. E^q 2 ls 
the energy per TiC>2 unit in the bulk oxide, E^f lk is the energy per Ti atom in the bulk metal, and 
2S is the area of the two interfaces in the interface structure. 

The results are summarized in [8l The negative values indicate that the interface is favored 
relative to the bulk. In contrast to the results for the O-O-Ti terminated interface shown in [5J 
' HCP-2' stacking here is less stable than the ' TT' stacking by about 0.1 J/m 2 . 

To illustrate the possible effect of the free surface on the W a d obtained by the film-like model, 
in ITOl we display the structure of the Ti-O-0 free surface from the film-like interface slab with 
' TT' and ' HCP-2' stacking, compared with its structure in the pure TiO2(100) calculation. After 
relaxation in the surface slab (flOf b)). Oi and O2 move significantly outwards, and finally, the Ti-O- 
O termination transforms to an arrangement like O-Ti-O. In Ti02, the surface terminated with Ti is 
highly polar and this outward movement of oxygen can decrease the dipolar moment. Very similar 
relaxation is observed in the interface with ' HCP-2' stacking (PTOld)). However, for the ' TT' 
stacking (flOTc)). the relaxation of Oi and O2, does not occur, and the termination remains Ti-O-O. 
As we mentioned above, in ' TT' stacking, the first Ti layer of the metal sits in a position like 
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the extension of the oxide, so the O-O-Ti terminated interface with a stoichiometric oxide can also 
be seen as a Ti-O-0 terminated interface with a non-stoichiometric oxide. This different surface 
relaxation directly affects the work of adhesion of the interface. We do not think this kind of 
interaction between the interface and polarized surface can be removed by increasing the thickness 
of the oxide slab. 
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Figure 10: Structure of the Ti-O-0 terminated TiO2(100) surface, (a) unrelaxed structure (bulk 
termination), (b) relaxed structure in the TiO2(100) surface slab, (c) relaxed structure in the O-O- 
Ti terminated interface slab with ' TT' stacking, (d) relaxed structure in the O-O-Ti terminated 
interface slab with ' HCP-2' stacking. Red spheres are O atoms, and silver spheres are Ti atoms. 



Table 7: Bader charge (e) of the selected atoms on Ti-O-O surface included in different 
structures, the labels of the atoms are identified in[10tb), (c), and (d). Bulk values for Bader 
charge 2.18<? and -l.lOe for Ti and O in bulk rutile. 

surface Ti02 unit inside Ti02 unit 



Atom 


Tii 


Oi 


2 


Ti 2 


3 


4 


free surface 


1.94 


-1.03 


-1.19 


1.88 


-1.23 


-1.13 


TT stacking 


1.56 


-1.20 


-1.16 


2.11 


-1.14 


-1.10 


HCP-2 stacking 


1.93 


-1.04 


-1.19 


1.88 


-1.25 


-1.14 



In order to analyze the effects of interfaces on the Ti-O-0 surface in detail, we examine the 
Bader charge on the ions near the surface in P7J In all cases the electron density near the surface 
is enhanced by almost one electron per surface cell. The charge distribution near the free surface 
is unaffected by an ' HCP-2' interface, but the presence of a ' TT' interface, changes the surface 
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charge distribution significantly. This interaction between surface and interface cast doubt on the 
validity of the film-like model. Similar calculations for the single-oxygen structures (terminated 
with O-Ti-O), show no such structure or Bader charge variation. Thus for the double-oxygen 
interface, the film-like model is still reliable. 

So get an estimate of finite size error, we can directly compare the calculated work of adhesion 
for the same interface obtained by bulk-like and film-like calculation d5]),and this requires subtrac- 
tion of the energy of the free surfaces. These are: 1.99 J/m 2 and 2.20 J/m 2 for Ti(1010) terminated 
with small and large interlayer spacing respectively; 0.68 J/m 2 for TiO2(100) terminated with O- 
Ti-O, and 7.91 J/m 2 for the sum of Ti-O-0 and O-O-Ti terminations. The resultant interface energy 
WfUm is listed in[8l For the single-oxygen interfaces, V^buik is m reasonable agreement with Wfu m . 
For the double-oxygen interface with ' TT' stacking the apparent discrepancy is obviously due to 
the different surface reconstruction, as just described. And the bulk-like model is more reasonable. 
In the next subsections, we will discuss the geometry as well as the electronic structure difference 
for these two models. 

Table 8: Relaxed values for W^k (J/m 2 ), compared with Wfu m that obtained by excluding the 
contribution of the surface energy from the work of adhesion by the film-like model. 



Termination 




single-oxygen 




double- 


■oxygen 


Stacking 


OT 


HCP FCC 


FCC-2 


TT 


HCP-2 


Wbulk 


0.30 


-0.03 -0.95 


-0.93 


-0.60 


-0.51 


WfUm 


0.39 


0.10 -0.82 


-0.88 


-0.05 





Atomic, electronic structure of single-oxygen interface 

The structure of ' FCC (fTTj) is quite similar as the ' FCC-2' ©, in particularly the configuration 
near the interface and relaxations of the interfacial atoms (Tii , Ti3, Oi , O2) in the relaxed structure. 
Atom spacing Oi-Ti2 and 02-Tii is 2.10 A and 2.12 A, very close to the results in 'FCC-2' 
stacking obtained by the film-like model, the Bader charge © and projected DOS on the interfacial 
atoms (fT2l) are also similar to the film-like model. For the total DOS @d)), the surface-atom- 
induced small peak around -17.5 eV is absent, which is expected. 
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Figure 11: Structure of single-oxygen interface with 'FCC stacking, viewed along [001] direc- 
tion of Ti02- Only layers near one interface are shown: the other is symmetrically equivalent. 
Left: without atomic relaxation (the equivalent atoms on the boundary are also shown). Right: 
with atomic relaxation. Red spheres are O atoms, and silver spheres are Ti atoms. 
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Figure 12: DOS projected onto selected atoms for the single-oxygen interface. The label of the 
atom is identified in[TTl 
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Atomic, electronic structure of double-oxygen interface using symmetric structure 
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Figure 13: Structure of double-oxygen interface with ' TT' stacking (due to the symmetry, only 
layers near one interface is shown), viewed along [001] direction of Ti02- Left: without atomic 
relaxation (the equivalent atoms on the boundary are also shown). Right: with atomic relaxation. 
Red spheres are O atoms, and silver spheres are Ti atoms. 



Not surprisingly. [131 is very similar to the ' TT' stacking of Ti-O-0 terminated interface ob- 
tained using the film-like model. For example, the layer spacing of Tii-Ti2 is 0.54 A, corresponding 
to 0.55 A for Tip Tis layer in [8l while the spacing of Ti2-Ti3 is 1.90 A, very close to the Tii-Ti2 
spacing 1.88 A ©. Note that Tii in[T3]is in the equivalent position to Tis m El 

Again the Bader charge © and projected DOS (fl4l) agree well with the film-like model ©. 
Comparing the total density of states Qc) and Oe)), only the surface oxygen states in the range 
from -17.5 eV to -15 eV are absent. 
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Figure 14: DOS projected onto selected atoms for the double-oxygen interface using symmetric 
structure. The label of the atom is identified in[ 
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Atomic, electronic structure of double-oxygen interface using asymmetric structure 
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Figure 15: Structure of two double-oxygen interfaces with different stacking sequence, viewed 
along [001] direction of Ti02- The ions in the center of the oxide are omitted. Left: without 
atomic relaxation (the equivalent atoms on the boundary are also shown). Right: with atomic 
relaxation. Red spheres are O atoms, and silver spheres are Ti atoms. 



As shown in \T5\ we can consistently reproduce the interface structure with ' TT' stacking 
([8] and [T3l) and ' HCP-2 ' stacking © very well in film, symmetric -bulk and asymmetric-bulk 
calculation. The Bader charges ([6]) are equivalent to film-like model. The relevant projected DOS 
compare very well with each other (fT6lvs[T4l. and[F7]vs|4]). Again no surface-atom-induced peak 
is found in the total DOS figure @f)). We can therefore be highly confident that these structures 
are converged with cell size and not artifacts of the boundary conditions. 



CONCLUSION 

First-principles calculation is performed to investigate the atomic structure and bonding nature of 
the Ti(1010) // TiO2(100)interface. A series of tests have been done to verify the settings used in 
the calculation, by studying the bulk and free surface properties and successful comparison with 
previous calculation and experimental work. 

Two different models are used to simulate the interface. With the film-like model, we consid- 
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Figure 16: DOS projected onto selected atoms for the ' TT' stacking interface. The label of the 
atom is identified in[T5l 
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Figure 17: DOS projected onto selected atoms for the ' HCP-2' stacking interface. The label of 
the atom is identified in[T5l 
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ered three terminations and different stacking sequences. Interfaces compatible with the extension 
mechanism (' TP stacking) were constructed for the O-O-Ti and Ti-O-0 systems, and we found 
these geometries quite favored, especially for the Ti-O-0 terminated interface. With the bulk-like 
model, different terminations and stacking sequences were also studied and the interface energy 
is calculated with no surface energy involved. Except in the one case with anomalous surface 
reconstruction, agreement between the two models is good. 

The relaxations tended to flatten the metallic and ionic surfaces. 

Substantial charge transfer in the interface structure is observed, giving interface charges ex- 
ceeding those in bulk Ti02 but this is confined to a small region around 5 A of the interface. As 
more electrons move to the oxide, the ionic bonding between the interfacial Ti atoms in the metal 
and the interfacial O in the oxide becomes stronger. These enhanced electrostatic fields at the inter- 
face are evidenced by the shift of the low-lying 2s state in the pDOS. More electrons are accepted 
by the interfacial oxygen ions than oxygens ions in the bulk rutile; this can be easily understood 
because the electron donators are more abundant in the presence of the metal slab. 

From the results of the projected DOS calculation, the picture emerges of metallic bonding 
changing to ionic bonding over a short region, with enhanced ionic bonding at the interface. In the 
O-O-Ti and O-Ti-0 interfaces, a small covalent component may exist, but for the Ti-O-0 interface, 
no hybridized states are observed on the interfacial Ti of the metal slab, thus no covalent interfacial 
bonding is formed. Also we notice that new gap states formed in the DOS of interfacial Ti of the 
oxide. As with structure and Bader charge, the atom-projected DOS is consistent between the two 
models. 

Most notably, the total energy of the slabs containing interfaces is lower than the equivalent 
numbers of atoms in pure Ti and pure Ti02- This shows that Ti has strong affinity for oxygen, but 
at the same time that the oxide is strongly bound to Ti. Combined with the good match of lattice 
parameters, this helps to explain why the oxide forms such a good protective coating for Ti alloys. 
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